Application of the multiple-scattering method to analysis of systems with semi-infinite 

photonic waveguides 
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We propose a technique of compensating the spurious reflections implied by the 
multiple-scattering (MS) method, commonly used for analyzing finite photonic 
crystal (PC) systems, to obtain exact values of characteristic parameters, such 
as reflection and transmission coefficients, of PC functional elements. Rather 
than a modification of the MS computational algorithm, our approach involves 
postprocessing of results obtained by the MS method. We derive analytical 
formulas for the fields excited in a finite system, taking explicitly into account 
the spurious reflections occurring at the artificial system boundaries. The in- 
trinsic parameters of the investigated functional element are found by fitting 
the results of MS simulations to those obtained from the formulas derived. De- 
vices linked with one and two semi-infinite waveguides are analyzed explicitly; 
possible extensions of the formalism to more complex circuits are discussed 
as well. The accuracy of the proposed method is tested in a number of sys- 
tems; the results of our calculations prove to be in good agreement with those 
obtained independently by other authors. 
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1. Introduction 

Photonic crystals (PCs) have recently become the ob- 
ject of increased interest as possible hosts for optical 
functional devices, e.g., beam splitters, demultiplexers, 
etc. [l[ These elements, as well as the basic build- 
ing blocks of photonic integrated circuits — waveguide 
bends, junctions etc. — have often been investigated by 
methods designed for finite systems, such as the finite- 
difference time-domain (FDTD) method [|J and the 
multiple-scattering (MS) method 0, 0|- Consequently, 
the devices in question were considered to be embedded 
in a finite fragment of a PC. This, however, involved spu- 
rious reflections at the artificial PC boundaries, signifi- 
cantly complicating the analysis of the device behavior. 

To remedy this situation, several methods suitable 
for analyzing systems with semi-infinite waveguides have 
been developed. In the effective discrete equations 
method [f| and the Wannier function method Q, the 
electric and magnetic fields are represented in a basis of 
states localized at elementary defects; various grating- 
based techniques, like those presented in @, H, 0|, have 
been developed as well. Both approaches, however, im- 
pose some restrictions on the problems to which they 
can be applied: the former is ill-suited to open sys- 
tems (i.e., those with vacuum regions extending to in- 
finity), while the latter only applies to systems with uni- 
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directional waveguides. Recently, the multiple multipole 
method has been extended by Moreno et al. [10] to sys- 
tems with semi-infinite waveguides; here, the fields on a 
transverse section of each waveguide, sufficiently distant 
from discontinuities (junction, bend etc.), are matched 
to a linear combination of the waveguide eigenstates. As 
demonstrated in [10j], this technique is very general and 
applicable also to the problems which cannot be dealt by 
the other above-mentioned methods. 

The aim of this paper is to show how to take into ac- 
count the presence of semi-infinite waveguides in calcula- 
tions performed on the basis of the MS technique. Owing 
to its particular simplicity and efficiency in dealing with 
the case of PCs composed of cylindrical rods, this tech- 
nique has gained significant popularity [Til fl2l lip . Al- 
though the approach proposed by Moreno et al. [1JJ] could 
be straightforwardly carried over to the combination of 
the MS method and the method of fictitious sources, for- 
mulated in this would require significant changes 
in the computational procedure, as well as availability 
of externally calculated data describing the fields corre- 
sponding to the waveguide eigenstates. In contrast, the 
approach outlined below consists solely in postprocess- 
ing of the results obtained by the 'pure' MS method. 
While admittedly less general than the other technique, 
it is nevertheless applicable to a number of situations 
frequently encountered in practice. Two of them are dis- 
cussed in Sections [5] and 02 in which the proposed proce- 
dure is applied to devices linked with one waveguide and 
two waveguides, respectively. Possible extensions of the 
method are discussed in Section [U 
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2. One waveguide 

A . Theory 

The first system to be considered is a single waveguide 
terminated at the surface of a semi- infinite PC [Fig .[Tfa)1- 
A frequently discussed problem (see, e.g., [l5L 1 1 6l. fl7^ ^) is 
the design of the precise shape of the waveguide outlet 
that would allow maximization of the power transmit- 
ted into free space. In this case, the main objective is to 
calculate the reflection coefficient of the outlet; the inten- 
sity of the field produced at some point of the free-space 
region (possibly at infinity) when an eigenmode of uni- 
tary power propagates towards the end of the waveguide 
is sometimes searched for, too. In the following we will 
show how both quantities can be found on the basis of 
calculations done for the finite system shown in Fig. QJb) , 
presenting a waveguide that is N unit cells long. 

The waveguide is assumed to be a single-mode one, and 
its unit cell to have a mirror-symmetry plane parallel to 
the yz plane, so that the propagating Bloch states of the 
waveguide are characterized by wave numbers k and —k. 
An arbitrary source S placed near the waveguide inlet ex- 
cites the right-propagating mode, which then undergoes 
multiple reflections. Let u n and d n (n = 0, 1, . . . , N) de- 
note the z component of the electric (TM polarization) 
or magnetic (TE polarization) field corresponding to the 
mode of wave number k and —k, respectively, at point n 
represented by coordinates (x = na,y = yo) in Fig. Hfb). 
The choice of j/o is arbitrary; the zero value can be as- 
sumed unless the waveguide eigenmode is antisymmetric 
with respect to the xz plane. From the Bloch theorem, 
we have 
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is the waveguide transfer matrix; a is the waveguide pe- 
riod and <P = e lka . The outgoing modes get partly re- 
flected at the waveguide ends; this can be expressed by 
the following 'boundary conditions': 



d N = ru N , 
u = u inc + r'd Q , 



(3a) 
(3b) 



where r and r' are the reflection coefficients at the wave- 
guide outlet and inlet, respectively, and u; nc stands for 
the effective field corresponding to the right-propagating 
mode excited by the source S, extrapolated to point 
n = 0. By combining Eqs. CD)-© and eliminating the 
variables uq, do, u n, and dpf, we get the linear system 
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Fig. 1. (a) The ideal system: a semi-infinite waveguide 
terminated at the surface of a PC. A single mode propa- 
gates along the waveguide from inside the PC; on reach- 
ing the waveguide outlet, it is partially reflected and par- 
tially transmitted into free space, (b) A finite counter- 
part of the system shown in (a), amenable to analysis by 
the MS method. The grey lines delimit the N unit cells 
of the waveguide. Smigaj-l.eps 



whose solution reads 



1 



J<p2N 



r $2N-r 



(5) 



In cells lying sufficiently far from the waveguide ends for 
the contribution of the evanescent states to be negligible, 
the total field f n can be expressed only in terms of the 
propagating modes: 



fn U n ~\~ d n 



2JV— n 



1 - rr ><p2N 



(0) 



We are interested in the reflection coefficient r, as well 
as m the intensity |/ ideal (<3)| 2 of the field which would 
be generated at some point Q in the free space if 
the waveguide were semi-infinite and the incident right- 
propagating mode carried unitary power. 
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The reflection coefficient can be found simply by least- 
squares fitting of the numerically calculated values of 
f n+i/fn to those calculated by the formula resulting from 
Eq. ©: 



71+ 1 



n+1 



2N-n-l 



fn 



r $2N-r 



(7) 



In addition to r, this gives also the value of <£, which can 
be used for calculating the wave number k. Let us note 
by the way that solving Eq. (T7]) for n = I — 1 and n = I, 
where I is some fixed integer, yields 



, , 1 fl-i + fl+i 

k = ± — arccos 

a 2f t 



(8a) 



and 



fi$ - fi-i 



$2(l-N)- 



(8b) 



these values, calculated for I ss N/2, i.e., near the wave- 
guide center, can be used as starting points in the nonlin- 
ear least-squares fitting procedure. The sign of k should 
correspond to the physics of the problem at hand. 

It should be stressed that only cells distant enough 
from the waveguide ends (i.e., those labeled n — B, B + 
1, . . . , N — B, where the 'margin' B is a sufficiently large 
integer) should be taken into account in the above fit- 
ting procedure, since Eq. ([7]) has been derived with the 
assumption that evanescent states are of negligible am- 
plitude in the cells labeled n and n + 1. 

Let us proceed to the determination of |/ ldeal (Q)| 2 - 
The intensity \f(Q)\ 2 of the field generated at point Q in 
the finite system can be written as 



|/(Q)| 2 = |t| 2 |r(Q)| 2 M 



l*l 2 KQ)| 2 | 

|1 _ rr '$2N\2 I 



(9) 



where the transmission coefficient |i| 2 = 1 — |r| 2 repre- 
sents the fraction of the total energy emitted into free 
space when a waveguide mode reaches the outlet, and 
the 'transfer coefficient' t(Q), dependent on the position 
of Q and the geometry of the waveguide outlet, but not 
on the waveguide length, is a measure of the amount of 
this energy getting to point Q. In the ideal case of the 
semi-infinite waveguide, r' would be zero; thus, 



|/ Ideal (Q)| 2 = |^|r(Q)| 2 |«^ al | 2 



(10) 



with |^inc aL | 2 chosen so that the incident mode would 
carry power i? ld c cal — 1. From Eqs. ([9])-([T0"]) we get 
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since the incident power in each case is proportional to 
Ui nc | 2 . As shown in the Appendix, the total power flow- 
ing through an arbitrary transverse section of the wave- 
guide is equal to 



P = 



1 



1 



(12) 



This power can be easily calculated by the MS method. 
By using Eq. (|12[) for eliminating from Eq. I|ll|) the factor 
1 — rr'<P 2N \ 2 , which contains the unknown coefficient r' , 
we finally get 



|/ ideal (Q)l 2 = (i- 



pidcal 

>^H/(Q)I 2 - 



(13) 



To sum up, the data obtained in a single MS calcu- 
lation performed for a finite system with an A^-cell-long 
waveguide suffice for determination of the reflection coef- 
ficient of the waveguide outlet, as well as the 'corrected' 
field intensity in free space. When only the former quan- 
tity is required, the field excited at A^+l sites lying along 
the waveguide axis is all that needs to be calculated in 
the simulation; otherwise, the power flow through an ar- 
bitrary transverse section of the finite waveguide must be 
computed too. 

B. Examples 

As an example of application of these results, let us first 
consider the system shown in Fig.[lja): a Wl-type wave- 
guide embedded in a hexagonal lattice of dielectric cylin- 
ders with permittivity 11.56 and radius 0.18a, where a is 
the lattice constant. The surface termination creates a 
slight tapering of the waveguide exit. Figure [5] presents 
the frequency dependence of the reflection coefficient |r| 2 
in this configuration; the results obtained by the pro- 
posed technique are compared to the data reported in [l8j 
[Fig. 7(c)], acquired by the plane- wave-based transfer- 
matrix method. Our calculations were done for wave- 
guide length N — 15 with margin B = 4. The con- 
vergence is very fast: in fact, the data obtained in the 
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Fig. 2. The reflection coefficient of the waveguide outlet 
shown in Fig. [TJ Circles: the calculation results of Lin 
and Li [18l|: line: results of our calculations. Smigaj-2.eps 



4 



N = 7, B = 2 case (essentially the shortest waveguide 
for which fitting makes sense) differ by at most 6% from 
those plotted in the graph, and for frequency values below 
0.468 x 2irc/a the difference does not exceed 1%. There 
is a good agreement between our results and those of Lin 
and Li [l8j |. except for a small peak at frequency value 
u> w 0.46 x 2nc/a, present in the |r(w)| 2 plot obtained 
by these Authors, but not reproduced by the curve re- 
sulting from our calculations. Since the waveguide is a 
single-mode one and the mode dispersion curve is smooth 
around this frequency value, the physical origin of this 
sharp peak is not clear to us and we believe it might be 
a numerical artifact. 

As a second example, let us consider the leaky- wave 
photonic antenna 0, HH shown in Fig. [SJa). Here, 
the surface surrounding the waveguide outlet is corru- 
gated with period 2a and supports leaky surface modes, 
which are excited by the radiation coming from the wave- 
guide. At a certain frequency value waves scattered 
at individual perturbed surface cylinders interfere con- 
structively along the surface normal to produce a colli- 
mated beam. Figure J3Jb) presents the frequency depen- 
dence of the far-field intensity of the radiation emitted 
perpendicularly to the surface; the outcomes of 'naive' 
calculations done for finite-length waveguides are jux- 
taposed with those obtained by the proposed scheme 
for the semi-infinite waveguide. All curves are normal- 
ized to their absolute maxima. Evidently, in the finite- 
waveguide case, the shape of the main peak depends 
on the precise length of the waveguide; this dependence 
is especially pronounced to the right of the maximum, 
where the intensity values obtained for waveguides 11 
and 12 cells long differ by as much as a factor of two. 
On the other hand, the sidelobes in the frequency range 
uj < 0.40 x 2-Kcja are not significantly affected by changes 
in the waveguide length; this leads to the conclusion that 
they result from interference occurring at the crystal sur- 
face rather than in the waveguide. 

3. Two waveguides 

A . Theory 

In this Section, we will focus on the more complex system 
with two linked semi-infinite waveguides; the link can be 
realized, for instance, by a junction, a bend or a resonant 
cavity where scattering can occur. Our goal is to calcu- 
late the reflection and transmission coefficients of this 
discontinuity. The waveguides are again assumed to be 
single-mode ones and to possess a transverse symmetry 
plane; to these assumptions let us add that of their iden- 
tical geometry. Figure 0] depicts the finite system used 
in numerical calculations, where the left and right wave- 
guide comprise N and M unit cells, respectively, and the 
field is excited by the source S. 

Let u n and d n denote the field (electric or magnetic, 
depending on the polarization), at point n = 0, 1, . . . , N 
inside the left waveguide in Fig. [51 corresponding to the 
incoming and outgoing guided modes, respectively. Their 
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Fig. 3. (a) A photonic crystal leaky-wave antenna. The 
radii of the bulk and the surface cylinders are 0.18a and 
0.09a, respectively, a being the lattice constant; all cylin- 
ders have permittivity 11.56. Every second surface cylin- 
der is shifted by 0.3a towards the bulk. The corrugated 
surface is n sur f = 9 unit cells long, (b) The far-field in- 
tensity of the radiation emitted perpendicularly to the 
surface, calculated for waveguide length N — 11, 12, and 
13, as well as for the semi-infinite waveguide on the basis 
of data obtained for N = 12, B = 3. Normalization: sec 
text. Smigaj-3.eps 

right- waveguide counterparts at point m = 0, 1, . . . , M 
will be labeled u m and d m . The junction linking the 
waveguides can be described by its scattering matrix S 
defined by 



d_N 


= S 


UN 


S = 


P T 


d M _ 








7 P. 



(14) 



p, p and r, f denote the respective reflection and trans- 
mission coefficients of the junction. The 'boundary con- 
ditions' are in this case 

uq = rdo, (15a) 
u = u inc + rd , (15b) 

where r is the reflection coefficient at the waveguide outer 
ends (the waveguides having the same geometry, their 
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terminations can be assumed to be identical too), and 
Ui nc is defined as in Section O The Bloch theorem gives 



n = 1 



N 



M 



— 2 
1 

m — 



Fig. 4. An example system with two waveguides linked 
by a junction (in this case, a bend). Grey lines delimit 
the unit cells of both waveguides. Smigaj-4.eps 
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(16) 
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with the waveguide transfer matrix T defined, as before, 
by Eq. By eliminating the variables uq, do, uq and do 
from Eqs. (fl"5|) and (fT6|) , and the variables un, dpj, um, 
and dM from Eqs. (fl~4"|) and (jTTJ) , we obtain the linear 
system 








the solution of which yields the expressions for the total fields in the left and the right waveguides: 
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(18) 



(19a) 
(19b) 



The above procedure can be repeated for the situation 
when source S is replaced with source S located near the 
entry of the right waveguide; in this case, the fields g n 
and g m in the left and the right waveguides are given by 



g n = 



f(r<P n + <p- n )<p N + M 



(1 - rp$ 2N ) - r[p + r(rf - pp)$ 2N ]<P 2M 



(20a) 



(1 - rp<P 2N )<P m + [p + r(rf - pp)<P 2N }& 



2Nlsf,2M-m 



(1 - rp<P 2 



[p + r(rf - pp)<P 2N ]<P 2M 



(20b) 



Mi nc being defined analogously to u mc . 

It is important to note that the MS method allows very 



efficient calculation of the fields excited in a single struc- 
ture by several independent sources (e.g., S and S), since 
the scattering matrix of the whole system, whose diago- 
nalization is by far the most time-consuming step of the 
computational algorithm, is independent of the incident 
field Jl, i]. In the following we show how the values of 
fn, fm, 9n, and g m , calculated by the MS method, can 
be used for determination of the junction parameters p, 
p, r, and f . As before, in all these computations only 
cells sufficiently distant from the waveguide ends should 
be taken into account. 



1. We have 



hn+l 



r ^m+l _|_ 



r <pr. 



(21) 
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therefore, as in the one-waveguide case, the param- 
eters & and r can be obtained by least-squares fit- 
ting of the above formula's right-hand side to sim- 
ulation results. Good starting points for the fitting 
procedure are in this case 



, , 1 fl-i + fi+i 
k = ± — arccos = 



(22a) 



of six equations with six unknowns: fx, fi, v, v, u- v 
and u\ nc . Its solution reads 



j" = 



OL20LJ, 



OtlOtH — OL^OL-i 

a 2 (a 4 — r/3 4 ) 
cticti — aids ' 
"inc = olx - r(3i, 



v 



LI : 



rai/34 



a 3 (ai - r[3i) 



(28a) 
(28b) 



OL\Ct4 — OI2OL3 

= (X4,- rp4. (28c) 



and 



fl-l$ - fl 



(22b) 



with I w M/2. 



2. To shorten the notation, we introduce the following 
symbols: 



p = rp$ 21 \ Ji = r0 IM , 

r2N -*2M 



(23) 

v = rr&" , v = rf <£ ZJW , (24) 

C = Lip, - vi>, r) = l-Li-p, + (; (25) 



consequently, the formulas for /„, f m , g n , and g ri 
become: 



_ (l-/i)<P n + r- 1 (/x-C)<Z ) -" 

In ^inc) 
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-$ N - M u w 



(1 - Li)$ m +r- 1 (/2-C)^" 
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(26a) 
(26b) 
(26c) 
(26d) 



Thus, linear least-squares fitting of the numeri- 
cally calculated values of /„ to the function ai<P n + 
f3\<I>~ n allows to find the coefficients 



1-/2 

Oil = u- u 

V 



01 



rj r 



(27a) 



Similarly, fitting the values of f m , g n , and g m to 
the functions (a 2 <£ m + 2 $~ m )$ M ~ N , (a 3 $ n + 
f3 3 <p- n )<P N - M , and a 4 <Z> m + 04$~ m , respectively, 
yields the values of the coefficients 



OL4 



3. With definitions (f2l>)) of £ and rj included, the for- 
mulas for ai, 0i, a 2l Q!3, o;4, and /?4 form a system 



4. The reflection and transmission coefficients p, p, 
r, and f can now be calculated from Eqs. (|23|) 
(|24|). since the values of all the other parameters 
are already known. 

B. Example 

To test the accuracy of our method, we are going to 
apply it to the waveguide bend depicted in Fig. [4] A 
number of numerical studies of this system, using dif- 
ferent methods, are available in the literature pH la. ITc| . 
providing a reference for new techniques. The super- 
imposed plot in Fig. [5] shows the frequency dependence 
of the reflection coefficient |/o(cj)| 2 (since the bend is 
symmetric, p = p) calculated according to the above- 
presented scheme against the spectra obtained by the ef- 
fective discrete equations method [5|, the Wannier func- 
tion method [|| , and the multiple multipole method [l(| • 
Clearly, the overall shape of all four curves is very simi- 
lar, although the exact values of |p(a>)| 2 differ, especially 
near the boundaries of the bulk crystal gap. The spec- 
trum resulting from our calculations is almost identical 
to that calculated by Moreno et al. [l(| This may be 
partly due to the similarity of the MS and multiple mul- 
tipole methods (in both of them the field is expanded in 
a basis of Bessel and Hankel functions); however, it must 
also be noted that these methods are intrinsically 'exact' 
in the sense that the Maxwell equations are solved rig- 
orously, which makes the calculations accuracy only de- 
pend on the maximum order of the basis functions kept, 
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5. The reflection coefficient of the bend shown in 
Fig. |4] calculated by different methods (see text). The 
data plotted with circles, crosses, and squares have been 
taken from @ and respectively. Smigaj-5.eps 
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and the results are known to be rapidly convergent with 
the truncation order 4]. In contrast, the derivation of 
the effective discrete equations [f| involves a number of 
approximations (for instance, only the monopole eigen- 
modes of the elementary defects are taken into account). 
This may be the reason why the reflection coefficient cal- 
culated by Mingaleev et al. takes values ~ 30% larger 
than those obtained in our approach. 

4. Extensions 

The two device classes discussed above, i.e., those linked 
with one or two semi-infinite waveguides, represent a 
large fraction of commonly investigated photonic build- 
ing blocks. Though the two-waveguide case has been ana- 
lyzed with the assumption of identical geometry of both 
waveguides, this condition can be eliminated, with the 
transmission coefficients renormalized to refer to eigen- 
modes carrying unitary power, in a way similar to that 
presented for a single waveguide at the end of Sect ion [2X1 
Our formalism can be extended to systems with more 
than two waveguides as well. However, many such sys- 
tems, including most T and Y junctions discussed in the 
literature, are symmetric with respect to the axis of one 
or more constituent waveguides. Therefore, they can also 
be analyzed with the simpler two-waveguide formalism 
presented in this paper. To see this, note that the scat- 
tering matrix S of a three-waveguide junction is defined 
by [cf. Eq. mm] 







U N 




= S 











with S 



Pi 721 T31 
Tl2 P2 T32 
Tl3 T23 P3 



(29) 



where u l N and d l N refer to the incoming and outgoing 
mode, respectively, in the ith waveguide (i — 1,2,3). 
For simplicity reasons, all the waveguides are assumed 
to be of the same length N. If waveguides 2 and 3 are 
identical and the axis of waveguide 1 is a symmetry plane 
of the system, we have p 2 = p3, T12 = i~ 13 , r 2 i = T31, 
and T23 = T32. To these equalities we can add u 2 N = 
and d 2 N — dj^, provided that the sources are placed 
symmetrically with respect to that axis. As a result, 
Eq. ((291) simplifies to 



"4" 




Pi 2r 2 i 






T12 P2+T 23 



(30) 



Therefore, the junction can be treated as a two- 
waveguide system with effective reflection and transmis- 
sion coefficients given by the above formula. The val- 
ues of these parameters can be calculated by the method 
presented in Section [3l Evidently, the coefficients p 2 and 
T23 occur only in the form of their sum, and therefore 
cannot be obtained independently However, normally, 
the quantities of interest are the coefficients p\ and ri 2 , 
which are related to transmission and reflection of the 
radiation incoming from the first ('input') waveguide. 
Equation (|3"0"1) clearly shows that these parameters can be 



obtained straightforwardly from the two-waveguide for- 
malism. 

5. Conclusions 

We have presented a method, based on multiple- 
scattering numerical simulations performed for finite sys- 
tems, allowing to find the reflection and transmission co- 
efficients of photonic crystal functional elements linked 
with ideal semi-infinite waveguides. As our approach 
does not involve modification of the existing code imple- 
menting the multiple-scattering method, no serious pro- 
gramming and testing effort are necessary for its appli- 
cation. The proposed formalism allows for dealing with 
a wide variety of photonic crystal building blocks likely 
to be encountered in practice, as demonstrated by its 
successful application to a tapered waveguide outlet, a 
photonic crystal leaky-wave antenna, and a waveguide 
bend. 

Appendix: Power flow in finite waveguides 

This Appendix presents the derivation of formula (j 1 2 j) for 
the power flowing through an arbitrary transverse section 
of the waveguide shown in Fig. [IJb) . Only TM polariza- 
tion will be considered; the procedure to be followed in 
the case of TE polarization is completely analogous. 

The power P(na) flowing to the right through a plane 
x = na (perpendicular to the waveguide axis) per unit 
length in the z direction is given by the integral of the 
x component of the time-averaged Poynting vector S: 



P{na) 



S x dy 



E z (na,y) H*(na,y) dy 



(A.I) 



From Eq. ([(jj) we have 



E z (na,y) 



2N-n 



1 - rr ><p2N 



-e z (y), 



(A.2) 



where e z (y) denotes the electric field corresponding to 
the right-propagating waveguide eigenmode of amplitude 
satisfying the condition e z (0) = Ujn C . Similarly, for the 
magnetic field we obtain 



H y {na,y) = 



2N-r 



1 



t$2N 



-h y (y); 



(A.3) 



note the change of sign in the numerator, which results 
from the magnetic field being a pseudovector. By includ- 
ing these formulas into Eq. (|A.1|) . after some straightfor- 
ward algebra, we arrive at 



P{na) 



1 



|l-r 
x{(l 



! Re J - 2 Im[r<£ 2(A, - n) ] Im j} , 

(A.4) 



8 



where 



e z (y) K{y) dy. 



From the conservation of energy, P(na) must be constant 
for all n, since the waveguide is not leaky. We conclude 
that the second term in braces in Eq. (|A.4[) , being n- 
dependent, must be identically zero; this leads to the 
following condition for the fields corresponding to the 
waveguide eigenmode: 



Im J = 0. 

In consequence, Eq. (|A.4|) simplifies to 
\l-rr'$ 2N \ 2 nc ' 



(A.6) 



(A.7) 



where Pi nc = Re J denotes the power carried by the in- 
cident mode. 
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